Effect of microwave radiation on non-linear resistivity of a two-dimensional electron 

gas at large filling factors 
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We study transport properties of a two-dimensional electron gas, placed in a classically strong 
perpendicular magnetic field and in constant and oscillating in-plane electric fields. The analysis 
is based on a quantum Boltzmann equation derived for a weakly disordered two-dimensional elec- 
tron gas. We consider disordered potential with both long and short range correlations. Electron 
scattering off such disorder is not limited to small change in momentum direction, but occurs on 
an arbitrary angle, including the backscattering. The non-linearity of the transport in the consid- 
ered system is a consequence of two co-existing effects: formation of a non-equilibrium distribution 
function of electrons and modification of the scattering rate off the disorder in the presence of dc 
and ac electric fields. This work describes both effects in a unified way. The calculated dissipative 
component of electric current oscillates as a function of the electric field strength and frequency of 
microwave radiation in a qualitative agreement with experiments. 

PACS numbers: 73.23.-b, 73.40.-c, 73.50.Fq 



I. INTRODUCTION 

The discovery of the microwave-induced resistance 
oscillations (MIRO)^'^ and the zero resistance states 
(ZRS)'^'^ has raised interest in non-linear transport prop- 
erties of two-dimensional systems (2DES) in a perpendic- 
ular magnetic field at large filling factors. The dissipa- 
tive dc magneto-resistance exhibits giant oscillations with 
the inverse magnetic field when exposed to a microwave 
radiationi^i^i^ili^i^iia The period of MIRO is controlled 
by the ratio of the microwave frequency uj to the electron 
cyclotron frequency uJc = \e\B/mc in magnetic field B. 
In the high mobility samples the MIRO evolve into the 
ZRS when the linear dc resistance becomes negativejii 
These remarkable findings made it imperative to under- 
stand the mechanism of oscillations preceding the onset 
of the ZRS. The experiments were performed at rela- 
tively high temperatures and low magnetic fields, when 
the Shubnikov-de Haas oscillations arc suppressed. 

The appearance of the MIRO has been first attributed 
to the modification of the impurity scattering rates in 
the presence of the magnetic field J^ii^ii^ This scenario is 
commonly referred to as the "displacement" mechanism. 
It has been predicte d^^i^^ that a different, so-called "in- 
elastic" mechanism, dominates in the regime, where both 
MIRO and ZRS were observedi^ According to Ref. [H, 
the microwave radiation is responsible for formation of a 
non-equihbrium component of the distribution function, 
which is isotropic in momentum and oscillates as a func- 
tion of energy. The amplitude of such non-equilibrium 
component of the distribution function is characterized 
by the temperature dependent rate of inelastic scattering 
processes I/tqc due to the electron-electron interaction. 
The analysis of Ref. [HI suggests that in weak electric 
fields and at sufficiently low temperatures the "inelastic" 



contribution from the non-equilibrium component of the 
distribution function to the linear-response dc resistivity 
is larger than the "displacement" contribution. 

A different series of experiments focused on measure- 
ments of the non-linear differential resistance in the ab- 
sence of microwave excitationi ^^i^^i^^i^°i^^i^^ In these ex- 
periments, the differential resistance has been measured 
in the Hall bar geometry as a function of the applied di- 
rect current. This current creates a strong electric field 
in a perpendicular direction, known as the Hall field, 
provided that the magnetic field is strong, lUcTh- 3> 1, 
where rtr is the transport scattering time. The scat- 
tering off disorder in the Hall field is accompanied by 
a change of electron kinetic energy and leads to depen- 
dence of transport characteristics on the strength of this 
field. In particular, the differential resistance exhibits os- 
cillations, called the Hall induced resistance oscillations 
(HIRO), as a function of the Hall electric field strength 
E. The HIRO were explainedii as a result of the geo- 
metric resonance in the electron transitions between the 
tilted Landau levels when the diameter of the cyclotron 
trajectory becomes commensurable with the spatial mod- 
ulation of the density of states. More rigorous approach^'^ 
employing the quantum kinetic equation showed that the 
"inelastic" mechanism is important in a relatively narrow 
interval of applied electric fields and the "displacement" 
mechanism becomes dominant in the regime of strong 
direct current, where HIRO were observed. 

The effect of the microwave irradiation on the non- 
linear transport was experimentally investigated in 
Refs. J.4.25, in which a 2DES was subject to both con- 
stant and oscillating electric fields. The value of the dif- 
ferential magneto-resistance depends on two dimension- 
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less parameters 



where E is the magnitude of the in-plane constant electric 
field, i?c = vf/uc is the cyclotron radius, vp is the Fermi 
velocity and tUc is the cyclotron frequency; throughout 
this paper we use h = 1. Maxima of the magneto- 
resistance in the vicinity of the main diagonal of the two 
dimensional (eac, edc)-plane are obtained when the sum 
Eac + Cdc is integer. Interestingly, this simple rule does not 
hold farther away from the main diagonal e^c Cac- And 
in fact, the interplay between both types of excitation 
gives rise to an unexpectedly rich structure of extremes 
and saddle points of the differential magneto-resistance 
in the (cac, edc)-planei^ 

Oscillations of the differential resistance as a func- 
tion of edc are understoodi^'^^iii^'^^ in terms of elec- 
tron backscattering off impurities, which corresponds to 
change of electron direction on its opposite. Therefore, 
an appropriate model of disordered potential has to in- 
clude processes of electron scattering on an arbitrary an- 
gle 6, including 9 = tt. Such potential has both long 
range correlations being responsible for small angle scat- 
tering, and short range correlations leading to backscat- 
tering. The proper treatment of the disorder potential 
with above properties requires a further extension of a 
kinetic theory of 2DES,"'^"' developed for smooth disorder. 

The goal of the present paper is to construct a sys- 
tematic theory of magneto-oscillations of the differential 
resistivity of the 2DES in the presence of electric fields 
of arbitrary strength. In experiments of Refs. [23l25l . 
oscillations have been observed at large filling factors, 
~ E-p/u>c S> 1 with E-p being the Fermi energy. In 
this limit, we can treat the kinetics of electron gas semi- 
classically. The analysis is performed in the experimen- 
tally relevant range of classically strong magnetic fields, 
l/^tr ^ Wc- In the present work we focus on the situation 
when Landau levels are not resolved, which implies the 
inequality Wc^q ^ 1, where Tq is the quantum scattering 
time. The temperature is assumed to be relatively large, 
T > Wc, so that the Shubnikov-de Haas oscillations are 
exponentially suppressed; fee = 1. 

The two most important contributions to the non- 
linear electric current are the "inelastic" contribution 
originating from the modification of the electron distri- 
bution function^ii^i^'^ and the "displacement" contri- 
bution representing the changes of electron scattering 
amplitudes off disorder. Additional contributions were 
identified and studied in Ref. 26, but these contributions 
have additional smallness in systems with mixed disorder. 
The displacement mechanism can be studied by various 
methods,— ii^iii^i^liSS'^g -^^j^jc]^ provide a qualitatively cor- 
rect picture for electron transport in strong electric fields. 
However, for a full description of the crossover from weak- 
to strong- fields, the kinetic equation is necessary. 

The paper is organized as follows. In Sec.[lTlwe present 
a simplified analysis of magneto-oscillations in combined 



constant and oscillating electric fields and summarize the 
main results. In Sec. IIIII the kinetic equation is derived 
in the framework of the Keldysh formalism. We solve 
the kinetic equation within a bilinear response in mi- 
crowave field and apply this solution to calculation of 
the non-linear current in Sec. IIVI Sec. |V] contains an 
analysis of the current beyond the bilinear in microwave 
field response. Discussion and conclusions are presented 
in Sec. ED 



II. MAIN RESULTS 

A. Bilinear response to the applied microwave 
radiation in strong dc electric field 

In this section we present heuristic discussion of the 
results of the paper for the dissipative current in strong 
dc electric field, but consider the contribution to the elec- 
tric current that is only bilinear in the microwave electric 
field. Our analysis employs the semiclassical treatment 
of electron motion in crossed electric and magnetic fields, 
valid in the limit of high Landau level index E-p/ujc 3> 1, 
where Ep is the Fermi energy and ujc is a cyclotron pe- 
riod. According to this picture electron scattering off 
impurities amounts to the spatial shift of the guiding cen- 
ter of the cyclotron orbit. Fig. [1] In the presence of the 
electric, E = e^E, and magnetic, B = e^B, fields, the 
dissipative current results from the imbalance between 
the drift of cyclotron orbits parallel to the electric field. 
We denote a unit vector forming angle ip with the direc- 
tion of the electric field by n^p = {cosi^; s'mip; 0}. The 
electron scattering resulting in the change of the direc- 
tion of motion from to nip> leads to the shift of the 
guiding center AR given by 

AR^^^p> = RcBz X (n^ - n^>) . (2.1) 




FIG. 1: Scattering off impurity leads to the shift of the guiding 
center of a cyclotron electron trajectory. 

We present the current as a sum of two contributions 
J=ji+J2. (2.2) 
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Here, the first term 



X [/(e, x) - /(e, a; + e^l^R^^^,)] 



(2.3) 



describes ttie current in the absence of processes changing 
electron energy by absorption or emission of a microwave 
field quantum with energy uj. This term contains virtual 
processes of electron scattering in microwave field, which 
modify momentum scattering rate off disorder. Thus, 
function V]^^^, is the disorder scattering rate for the di- 
rection of electron momentum from n,^ to n^/ and has 
the following form within the Born approximation: 



1 



i^{e,x + AX^^^,) 



^0 



(2.4) 



where is the disorder scattering rate in the ab- 

sence of electric and magnetic fields, and 



1 



1 — cos((^' — ip) 



The dimensionless parameter 



,,2 „2 p2 



(2.5) 



(2.6) 



is proportional to the microwave power. The upper 
(lower) sign in Eq. (|2.6p corresponds to the right (left) cir- 
cular polarization of the incoming microwave radiation. 
Close to the cyclotron resonance the dynamical screening 
of the microwave field by electron system has to be taken 
into account. This screening results in a modified form 
of parameter y^^, see Ref. |30- The drift of the guiding 
centers due to the microwave field effectively smears the 
disorder potential felt by the electron. For that reason 
the scattering rate F^^^,, Eq. (|2.4p . is suppressed by the 
microwave radiation. This scattering rate suppression is 
reminiscent of the "motion narrowing" phenomenon and 
is further discussed in Sec. IIII CI 

Apart from the virtual processes of electron interac- 
tion with microwave field, there are real processes, which 
are accompanied by the absorption and emission of pho- 
tons. The second term in Eq. (|2.2p takes into account the 
contribution to the electric current from such processes: 



J2 



± 



X [f{e,x)- f{e±Lo,x + e,AR^^^,)] . (2.7) 

The rate of such scattering processes can be written in 
the form 



F 



± 



iy{e±uj,x + AX, 



(2.8) 



We emphasize that in the present subsection we do not 
consider multi-photon processes. This approximation is 



justified if the dimensionless parameter CP^^ <C 1. We 
analyze the case of the arbitrary parameter in Sec.lVl 
see also Sec. IIIBl 

The angular integrals in Eqs. (|2.3p and (|2.7p are re- 
stricted by the conditions AX^^^pi — e^AR^^^i > 
0, and spatial integration is limited to the stripe 
—AXip^ip' < X < 0. The electron density of states has 
spacial dependence in constant electric fields. This de- 
pendence appears under the spacial integral and we dis- 
cuss this dependence in more detail. In a perpendicular 
magnetic field and in the absence of electric fields, the 
density of states of a 2DES, i^(e), has an energy modula- 
tion with period lUc 



i^{e) = i/Q 1 — 2A cos 



2TTe 



(2.9a) 



where i^q is the density of states in the absence of fields, 
and the factor A = exp(— Tr/wc^'q) ^ 1- A constant elec- 
tric field tilts electron density of states along its direc- 
tion, , resulting in spacial dependence of the density of 
states: 



iy{s, x) = v[e + eEx) . 



(2.9b) 



We present the results for the non-linear current in 
strong dc electric fields, edc = 2 |e| ERc/uOc ^ 1, where 
the differential resistance exhibits an oscillatory behav- 
ior. Rigorous analysis of Sec. IIV CI shows that in this 
case the dominant contribution to the non-linear current 
originates from the smooth component of the distribution 
function, which can be taken as the Fermi-Dirac distri- 
bution function at temperature T: 



Ms) 



1 



e^/T + 1 ■ 



(2.10) 



Substitution of Eq. ([^^ to Eq. ([^^ gives the contri- 
bution ji to the current in the form . 



Ji 



''■J'^Jde.ieHe^eEAX,^,,) 



1 



^ I, 



(2^) 
X AX,^ 

X [Me] ~ Ms + eEAX^^^,)] 



From Eq. (|2.1ip we obtain 



(2) 



(2.11) 



(2.12) 



Here the first term is the zeroth order in A, neglecting 
the oscillations of electron density of states in magnetic 
fields in Eq. (|2.1ip . This term corresponds to the classical 
Drude contribution to the current at large Hall angles, 



CTD 



2 Rc 1 
e 1^0 — , — 

TtT Ttr 



dOl- cos 6 



27r 



(2.13) 
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where rtr is the transport scattering time. The second 
term in Eq. (|2.12p represents a contribution to the cur- 
rent, which is non-hnear in the apphed electric field E 
and quadratic in parameter A: 



dipdip' 



(2^) 



CP 



[AX^p^^'] cos 



T(_p — (_p ' Tip -~ip' 

2'KeEAX,^^,^, 



(2.14) 



The contribution of the first order in A, which is omit- 
ted in Eq. p.l2p . describes the Shubnikov-de Haas os- 
cillations and can be estimated as cx ct^iAtA. Here the 
additional small prefactor At — exp(— 27r^r/a;c) appears 
as a result of averaging of the rapid oscillations in the 
density of states over thermal energies \e\ < T. The con- 
tribution of the second order in A contains a square of 
the oscillating component of the density of states and is 
not exponentially suppressed after integration over ther- 
mal energy window. In this paper we consider the limit 
of relatively high temperatures T ^ lOcI^-k. The latter 
condition is normally satisfied in experiments performed 
at the non-linear regime in dc electric field. Under this 
condition the quadratic in A contribution dominates over 
the Shubnikov-de Haas contribution linear in A. 

For Edc ^ 1 the angular integrations in Eq. (|2.14p 
can be performed in the stationary phase approximation. 
The main contribution comes from the scattering within 
narrow intervals centered at ip ~ ±'k/2, Lp' = 
These back scattering processes correspond to the 2i?c 
jumps of electron guiding center along the electric field. 
Equation (P7Ti|) yields^^ 



5jf^ « (2A)^ (1 - 2?^) 



sin 27redc 



(2.15) 



A similar analysis of the contribution Eq. (|2.7p to the 
electric current due to the real processes of absorption 
and emission of microwave photons gives: 



J2 



i{2\Yy, 



■ sin27r(edc + fac) 



■ sin27r(edc - Cac) 



(2.16) 



We emphasize that both Eqs. (|2.15p and (|2.16p are ex- 
pressed through the parameter edc oc ER^. These two ex- 
pressions are obtained within a stationary phase approxi- 
mation and correspond to processes representing shifts of 
electron cyclotron trajectories by distance 2Rc along the 
applied dc electric field. Therefore, Eqs. p.lSp and (|2.16p 
have a simple geometrical interpretation in terms of com- 
mensurability between the space modulation of electron 
density of states by electric field and the maximal dis- 
placement 2i?c in a single scattering process off disorder. 
The expression for the non-linear current beyond the sad- 
dle point approximation is evaluated in Sec. lIVBl 



In general, we represent the total current as a sum 
of the linear Drude term udE and the non-linear term 
which arises due to oscillatory density of states in 
perpendicular magnetic field: 



(2.17) 



For the overlapping Landau levels, when the density of 
states is given by Eq. (j2.9ap . the non- linear contribution 
is quadratic in parameter A. 

Within a bilinear response to the applied microwave 
electric field and in strong dc electric fields we combine 
Eqs. ([^ and ((TT^ and obtain 

5j-'5jf^+j2. (2.18) 

Equations (|2.15p . (|2.16p and (|2.18p represent some of 
the main results of this paper. As we discuss in 
Sec. llVi despite its simplicity, Eq. (|2.16p explains pro- 
nounced features of the transport measurement reported 
in Refs. 



Current at arbitrary microwave powers in 
strong dc electric field 



The result of Eq. (|2.18p is just a limiting case of ex- 
pression obtained in Sec. |V] for the current in strong dc 
electric fields. In this case we again can neglect the mod- 
ification of the electron distribution function by electric 
fields and consider only the displacement mechanism for 
generation of the non-linear component 5j of the current 
Eq. (|2.17p . We found the following expression 



5j .(2A)2»^ 



sin 27redc ^/o ( 4 V sin Treg 



cos 27redc cos Treac \/CP^ Ji f4\/5\^sin7rea 

Cdc V 



(2.19) 



Here J„(a;) are the Bessel functions. Performing an ex- 
pansion in Eq. (|2.19p to the first order in we recover 
Eq. (|TT5)1 in terms of Eqs. (|TT5|l and (PTB)) . 



C. Weak electric fields 

Above results were obtained under assumption that the 
distribution function of electrons is given by the Fermi 
distribution function, see Eq. (|2.10p . In the presence of 
electric fields the distribution function deviates from the 
equilibrium configuration. This non-equilibrium distri- 
bution function affects significantly the electric current 
in sufficiently weak electric fields. For smooth disorder, 
we recover the result of Ref. [l5| for the non-linear con- 
tribution bj to current in terms of dimensionless param- 
eters Edc, Cac, defined by Eq. (jl.ip . and defined by 
Eq. (IMD: 



8j = 2AVd£^- 



STreacCPt^ sin 27rea 



Ttr/7"ee + 7r2e^^/2-|-2T^sin TTEa 



(2.20) 
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Here I/t^c is the relaxation rate of the non-equihbrium 
component of the electron distribution function due to 
electron-electron interaction. In the presence of sharp 
disorder, the "displacement" contribution, arising due to 
the modification of electron scattering rate off disorder 
may become comparable to the "inelastic" contribution 
Eq. (|2.20p . The latter contribution survives only in rel- 
atively weak electric fields, edc 1- In stronger elec- 
tric fields, edc S> 1, it contains an extra small factor 
Tq/(r7redc) a-iid the "displacement" mechanism becomes 
more important, see Sec. IIVCI 



III. QUANTUM KINETIC EQUATION 

In the present section, we study a disordered 2DEG 
subject to in-plane electric fields and derive the quantum 
kinetic equation, following Ref. fl^ but consider a mixed 
disorder characterized by scattering amplitude which is 
finite for an arbitrary scattering angle. The Dyson equa- 
tion for the disorder averaged electron Green function 

(idt - G{tt') = lS{t -t')+ I dtit{tti)G{tit') , 



(3.1) 



where H is the one-electron Hamiltoniari^in the absence 
of disorder. Both the Green function G and the self- 
energy E for electron scattering off disorder are matrices 
in the Keldysh space: 



G = 





G^\ 






V 


G^ j 




I j 



(3.2) 



Matrix 1 stands for unit matrix in the Keldysh space as 
well as in the one-electron Hilbert space. 

Below we assume that the disorder potential is a 
combination of the long-range potential with correlation 
length ^ created by remote positively charged donors, 
and the short-range potential with much smaller correla- 
tion length. When conditions ^ ^ \h and p-p^ ^ 1 are 
satisfied, the self consistent Born approximation for the 
self-energy calculation is applicableiU, In this case, the 
expression for the self-energy takes the formM 



S(P) 



W{\q\) e'^'-G{p)e 



iqr 



(3.3) 



Here the function M^(|q|) is introduced as a Fourier trans- 
form of the correlation function of a Gaussian disorder 
potential 



(c/(n)c/(r-2)) 



W{q)e 



iq{ri—r2) 



(3.4) 



The external electric field can be eliminated by trans- 
ferring to the moving reference frame, r — > r -I- C(Oi 



where C,{t) describes a two-dimensional electron motion 
in crossed electric and magnetic fields: 



(3.5) 



Here ujc = |e| B/mc is the cyclotron frequency, and e is 
the antisymmetric tensor: eE — E x for any vector E 
lying in the x — y plane of the 2DES, see Fig. [TJ In the 
moving reference frame, the disorder potential acting on 
electrons is time dependent, and Eq. (|3.3[) becomes 



S(P) 



iqri 



where 



Wt,t,(g) = M/(g)e''^^'l'^ 



(3.6) 



(3.7a) 



(3.7b) 



To proceed further we introduce the operator of the 
guiding center coordinate R: 

iJ=|^-^, r = R+Xj,ep, (3.8) 
Zm 

where A^f — {ch/ \e\ BY^^ is the magnetic length. The 
commutation relations between the operators of the guid- 
ing center R and momentum p are: 

[Ra,Rp] = iXneap, \Pa,Pp] = -T2-£a/3, [R,P] = 0. 

(3.9) 

The equation (|3.6p takes the form 

S(p) = J ^W^2{qW'''P^"G{p)e-^'''P^"»]. (3.10) 

We first analyze the retarded and advanced compo- 
nents of the Dyson equation, Eq. (|3.ip , which determine 
electron spectrum. Then, we reduce the equation for the 
Keldysh component of Eq. (|3.ip to the kinetic equation 
for the electron distribution fmiction. 

A. Electron spectrum 



The Dyson equation, Eq. (|3.ip . for the retarded com- 
ponent of Green function is given by 



G^{t,h;p) 
S{t-h) 



27r 



dt2^'\t,t2)G''it2,ti;p) (3.11) 



along with the self consistency equation (|3.10p . We 
limit the analysis to the first order in the parameter 
A = e"'^/""'^''. The zeroth order solution in A corresponds 
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to the solution in the absence of a magnetic field. The 
standard answer for the self-energy in this case is 



2t, 



■IoS{t - t') , 



(3.12) 



where /g is the unity operator in the coordinate space. 

The Green function corresponding to (|3.12p is obtained 
by solving (|3.1ip 

G^{t, t') = ~ie{t - i')e-«-ff(t-t')e-(t-t')/2r, _ (3^^3) 

We consider the first iteration of the self consistent 
Born approximation. For this purpose, we substitute 
Eq. (|3.13p to the retarded matrix element of Eq. (|3.10p : 



(3.14) 



We notice that the operator product in the integrand in 
p.l4[) can be written as 



where the time dependent operators are pt = e"" "pe 



, (3.15) 



To find the time dependence of the operators, it is con- 
venient to introduce^ "raising," p+, and "lowering," p~, 
operators, defined as = {px ± ipy)/^/2 and obeying 
the commutation relation [p~^,p^] =^ The intro- 

duced notations for p^ are useful as new operators have 
simple form in the interaction picture 



(3.16) 



We also introduce the notations for vectors — [q^ ± 
i(7y)/V2 and write qep = q~ p^ — q'^ p~ . Using Eq. (|3.16p 
and the operator relation e^e^ 
further transform Eq. (|3.15p to 



gA+i3g(l/2)[A.B] 



= exp |2iA^((?f p+ + qjp ) sin 
X exp < sin [uct] 



OJct 


1 


2 





(3.17) 



q^e^"^"*^"^, or in 



In equation p.l7p . we denoted qf 
other words, the vector qt is obtained from the vector 
q by a rotation on angle ujct/2. Using Eqs. (|3.15p and 
()3.17p we can write Eq. (|3.14p in the following form 



d^q 



±'tf;{p) = -i6{t - t')e-{''ff+i/2r,)(t-t') 
X T^tt'(|q|) exp |2iA|^qt_t'Psin 
X exp Y sin [ujc{t - 



(27r)2 



(3.18) 



It is illustrative to consider the zero field limit. In this 
case the momentum operator becomes a c-number and 
Eq. ((3J8| yields 



t'i,{p)^-ie{t~t')e-^'-'"^'^^- 
d^q 



(27r) 



rW{q)e 



(3.19) 



Changing integration variables in p.l9p to ^p-q and the 
angle formed by the vector p — q with some fixed direc- 
tion, we show the consistency of the employed approxi- 
mation, namely 



with defined by Eq. ([3l^ and 



1 



de 1 

2tt Te 



The scattering rate 



1 



— = 2vq'kW 2pp sin - 

T0 \ 2 



(3.20) 



(3.21) 



(3.22) 



off disorder on angle 6 can be written in terms of its 
angular harmonics 



1 



-t-oo 

E 



n— — oo 



(3.23) 



Then, Tq = tq. 

In finite magnetic fields, we notice that the exponential 
factors in Eq. p.lSp are 27t/uJc periodic. The argument 
of these exponents vanishes at t — t' = ITc with integer 
I and the integral in Eq. (|3.18p diverges; Tc = 2tt/ll1c is 
the cyclotron period. We argue that this integral gives, 
in fact, a (5-peak of the width 5t ~ 1/E-p for time differ- 
ence t — t' — ITc- Indeed, this statement is obvious if the 
operators in the exponent of Eq. (|3.18p can be treated as 
commuting. For time intervals \t — t' — lTc\ < l/V-E-F^c 
the commutator of the two operators in the exponent 
of Eq. p.l8p is small since, in this case, each operator is 
multiplied by sva{L0c{t — t')/2) < yJuic/E-p and we can ap- 
ply the same argument as the one used in zero magnetic 
field. On the other hand, for 1//E^ <\t-t' - lTc\ < 
Tc/2 the result of the integration in ()3.18p vanishes be- 
cause of the rapid oscillations of the exponent. This can 
be checked explicitly by calculation of the matrix ele- 
ments of the self energy Eq. (|3.18p . We conclude that 
the non-commutativity of the operators can be ignored 
in Eq. (|3.18p and we can apply our zero field considera- 
tions whenever i — <' is a multiple of the cyclotorn period 

The contributions to the self-energy (|3.14p with I > 1 
are proportional to A' and, therefore, can be neglected 
in moderately weak magnetic fields, when A <C 1 . We 
stress that the terms with / > 1 are beyond the accuracy 
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of the first iteration of the self-consistent scheme. The 
corrections of the higher orders in A can be taken into 
account in the spirit of Ref.[3, where the small scattering 
angle was only considered. In this paper, we restrict our 
analysis to terms / = 0, 1, sufficient in not too strong 
magnetic fields. 

Introducing a new variable p' = p — q in (|3.18p and ne- 
glecting the variation of the function Wt^t' [p' — p) given 
by Eq. (|3.7p (see bellow) we can easily perform the in- 
tegration over the absolute value \p'\. This integration 
results in 



1 /I 



'-Sit-t')-X6it-t' -T,)]±t,t'M 



-'t,t';ip 



where we have defined the integral kernel 



(3.24) 



OCt,t';AF{v)} ^ I ^ 1 F{V') (3.25) 



2n 



with an arbitrary function F(ip)] = (cos 1^9; sin (p; 0). 
The negative sign of the second term in Eq. (|3.24p corre- 
sponds to the chemical potential ^ — sujc, with integer s, 
such that the density of states i^(/i) is at minimum. We 
note that the exact position of the chemical potential is 
of no importance for our final results. 

Equation p.24p is valid if the variation of the ma- 
trix element for the transitions in the moving reference 
frame, Wt,t' [p' — p) can be neglected in the course of in- 
tegration over the absolute value \p'\. This requirement 
leads to the limitation on the strength of electric fields. 
Since the relevant time scale is \t — t' — lTc\ ^ Tq <C Tc, 
which insures that the parameter A can be defined in 
Eq. p.24p unambiguously, the important range of the 
integration is \p'\ < 1/v-pTq. The requirement of smooth- 
ness of Wt^t' (p' — p) in the momentum variable means 
that the distance electron drifts over one cyclotron period 
must be smaller than the quantum length vpTq. In the 
case of a constant electric field, this condition amounts to 
vd/oJc ^ wpTq, or equivalently, edc ^ EpTq. If this con- 
dition is not met, the amplitude of oscillations as a func- 
tion of Cdc is reduced because of the finite broadening in 
time of the self-energy. In experiments)^'^ EpTq « 10^, 
Edc ^ 5 and the above condition is satisfied. However, the 
amplitude of oscillations in Ref. [H shows tendency to de- 
crease when Edc increases. In our model, this tendency 
can be accounted for by taking into account the depen- 
dence of the function Wt^ [p' — p) on the absolute value 
of momenta, which becomes stronger as edc increases. 



to the kinetic equation for the electron distribution func- 
tion /, related to the Keldysh component through 
the standard expression 



G 



K 



G' 



G^ - 2[G^f - /G-^] , 



(3.26) 



In the present analysis, the Wigner transformation 
f{t, t'; R,p) of the distribution function / in time and co- 
ordinate variables is independent of the "center of mass" 
coordinate R= {r + r')/2 due to the translational sym- 
metry. The peaked structure of the retarded and ad- 
vanced Green functions in Eq. (|3.26p makes f{t,e;R,p) 
to be independent of the absolute value of the momen- 
tum, p. The dependence on the direction of the mo- 
mentum is still to be retained. Notice that although the 
two components of the momentum operator are not com- 
muting, the momentum direction is well defined in the 
quasi-classical regime Ep 3> lOc. Therefore, we can write 
f{t,t';R,p) = f{t,t';nA = ft,t'-^p- The solution of the 
resulting quantum kinetic equation presented in Sec. IIVI 
is consistent with the assumptions made above. 

The distribution function fit, t'; n^) obeys the follow- 
ing kinetic equation: 



m + lH)-f] = Stdis{/} + Stec{/}, 



(3.27) 



where the notation [•; •] stands for the commutator in 
one-particle Hilbert space and the time variable. In 
Eq. (|3.27p the coUision integral Stcc{/} describes the 
electron-electron interaction and is discussed in the end 
of this section. The collision integral Stdis{/} represents 
scattering off disorder and can be written as the sum of 
scattering "out" and "in" terms 



Stdis — Stout 



Stin , 



where 



iStout{/} = P^'Z-Zs^], 

*Stin{/}- ^[S^'-S^ + S-^]. 



(3.28) 

(3.29a) 
(3.29b) 



Next, we analyze the collision integral in Eq. (I3.27p . 
We start our analysis with the scattering-" out" term, 
Eq. (|3.29ap . written as 



^Stouti/jt.t';!^ ~ j ^'td"-ipft" ,t'-<p — ft,t";ip^t",t'; 

(3.30) 

Using Eq. (j3.24p wc can perform integration over inter- 
mediate time t" in Eq. (|3.30p and obtain 



B. Electron distribution function 

In the previous subsection, we analyzed the effect of 
disorder on the spectral characteristics of electron Green 
function, determined by the retarded and advanced com- 
ponents of the Green function. Now we reduce the equa- 
tion for the Keldysh component of the Green function 



Stout{/}t,t';(p — —'^t,t;ip{^} ft,t':ip 
+ A (3Ct,t-Ta;ip{^}ft-T^,t'-:p + ft,t' -T^;:p^t' -T^,t' ■,:p{^}j , 

(3.31) 

where the operator 3Cf.t';<^{l} acts on a unity as defined 
by Eq. (jX^ and = iiT/ujc. 
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We now turn to the consideration of the scattering-" in" 
term, Eq. (|3.29bp . We express the self energies through 
the Green functions, using the self consistency condition, 
Eq. (|3.6p and the parametrization in Eq. (|3.10p 



X e' 



'iqepX 



« , (3.32) 



where the Green functions are given by Eq. (|3.13p . Com- 
muting the retarded (advanced) Green function with the 
left (right) exponent, wc rewrite (|3.32p as 



Stin{f}t,t' = i 



(27r)2 



dt"Wut'{q) 



Gf t„ e"i'P*-*" ^" ft",t' {p)e-"i'P^' 



't,t 



(3.33) 



Using e"^^P^^f/f,t'(p)e ^'J'^p^h = ft,t'ip^q), we commute 
the exponents with the distribution function: 



Stin{/}t 



_cfq_ 



dt"Wt.t'{q) 
"ft",t'{p-q) 



-ft,t"{p - q)e"'''P^»e-""P^'-*"^"G^, t, 



(3.34) 



Following the same line of arguments as in the derivation 
of ([^^ we put ((HTM)) into the form 



StiMKt',^ = / dt" 

X < I ^S{t - t") - X5{t - t" - T,) ) Xt,t'Aft".t';v} 
+ ( \5{t' - t") - \5{t' - t" - Te) ) Xt,t'Akt";v} 



(3.35) 



where "X-t^t'-.^fiftiMw} defined by Eq. (|3.25p . Perform- 
ing the time integration in Eq. (|3.35p we find 

— ^^t,t';ip{ft,t'-T^;ip} ■ (3.36) 



Collecting Eqs. p.3ip and (|3.36p and using i[H; ft,t';'^] = 
^cdipft,t';ip we can finally write down the kinetic equation 



d_ d_ 

dt ^ dt' 



d 



dip 



ft.t'-.ip — Stdis{/}t,t':ip + Stcc{f}t,t';ip 

(3.37) 

with the collision integral 

Stdis{/}t,t':ip = 'y^t.,t';ip{ft,t':ip} ~ '^^t.tiipi^} ft.,t':ip 

— ^'^t,t'[ip{ft-T^,t';ip} + ^'^t,t-T^;ip{^}ft-Tc,t'-ip 

— ^'^t.t' iipift.t' ~T^:ip} + ^ft,t' -T^-.ip'^t' -Ta.t';ip{^} ■ 

(3.38) 



The integral kernel Xt^t'-A^if)} defined by Eq. ((3?25)) 
is the generalization of the corresponding differential op- 
erator derived in Ref. for small angle scattering. 

We briefly discuss the term in the kinetic equation, 
Eq. p.37[) . representing the electron-electron interaction, 
Stee{/}- As it was shown in Ref. 15, electric fields pro- 
duce an isotropic non-equilibrium contribution to the dis- 
tribution function, which can only be stabilized by the 
inelastic relaxation mechanisms. In weak electric fields, 
the important inelastic relaxation is due to the electron- 
electron scattering. To take it into account, we keep the 
corresponding collision integral Eq. (|3.42p . which can be 
written in the energy representation for a steady in time 
distribution function /(e): 

Steo {/(£)}= / de' f dnM{n,e,e') 



J{e)f{e+)~f{e')f{e'_) - f{e)~f{e+)f{e')~f{e'_) 

(3.39) 

where /(e) = 1 - /(e), e+ = e + n, e'_ = e' - VL 
and M(ri, e, e') describes the dependence of the ma- 
trix element of the screened Coulomb interaction on the 
transferred energy O and the electron energies e and e'. 
The kernel M(r2, e, e') has been discussed in details in 
Ref. fill Below we use the linearized version of Eq. (|3.39p : 



Steo{/(e)} 



f{e) - Me) 



(3.40) 



Here, t^c is the inelastic relaxation time due to the 
electron-electron interaction, 

T2 



1 



oc 



47ret 



■In- 



max 



{T, ' 



(3.41) 



which can be obtained as a projection of the linearized 
electron-electron collision integral on the oscillating har- 
monic of the distribution function, see Eq. (|4.4ap below. 



C. Bilinear response in microwave field 

In this subsection we simplify the integral kernel %t.t' 
in the limit of weak microwave power, keeping only terms 
which are bilinear in the microwave electric field, i.e. lin- 
ear in power 7^, introduced in Eq. (12. 6|) . In the presence 
of microwave radiation 3Ct.t' has the oscillatory depen- 
dence on the time variable (i -I- t')/2. The distribution 
function in turn acquires non-stationary corrections os- 
cillating with the microwave frequency uj. It follows from 
the form of the kinetic equation, Eq. p.37p , that those 
corrections are small in the parameter of the order of 
l/cjTtr in systems with considered here mixed disorder. 
We therefore neglect their contribution to the distribu- 
tion function and consider only the stationary component 
of the distribution function. The latter can be found from 
the following equation 



i^cd^It-t'iv) = Stdis{/}t-t' + Stco{/}, 



(3.42) 
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where the symbol (...) stands for time averaging over 
one period of the microwave osciUations and we included 
the collision term due to the electron-electron interac- 
tion. Due to the 1-k jio periodicity of the collision integral, 
Eq. dSJll), the time -average of the kernel 3C, Eq. (j3.25p . 
is given by the integral over one period of the microwave 
field. 

We consider the response of a two-dimensional electron 
gas to the in-plane electric field 

-Etot = -E + £^inw , (3.43) 



Averaging Eq. (|3.48p with respect to time results in 



[( 



""•vOCtT.tJ =— 2-(ny-n^') sm . 

(3.49) 



And finally the collision kernel to the second order in the 
microwave field takes the form 



represented as a superposition of a constant electric field 
E and a circularly polarized microwave field 



(3.44) 



where we introduced the complex polarization vector 
— ± iey)/^/2 with the property ee^ = ±ie^. In 
the above equations, the upper (lower) sign corresponds 
to the right (left) polarization of the microwave field 
propagating in the magnetic field direction, see Fig. [T] 

The displacement Cti,t2i Eq. (j3.7b|) . is found by solving 
Eq. ((3?5|) with the electric field specified by Eqs. ((03)) 
and (|3.44p . The linearity of Eq. (|3.5p allows us to rep- 
resent its solution as the sum of the displacements in 
constant and microwave fields, 



where 



dc 



/■dc I ^ac, 



e{t2 - ti)eE 
mwc 



(3.45) 



(3.46) 



f l-cosuj{t-t') 



(3.50) 



where the rate 1/f has been introduced in Eq. (|2.5p and 
the quantity 



Wiptp' = eERc{smip — sinip') 



(3.51) 



is the work done by the dc electric field as the result 
of scattering off an impurity^^ W^^' = eEAR^^^i 
with the shift of the cyclotron orbit Ai?^^^/ given by 
Eq. ISTTli . 



The collision integral due to scattering off disorder is 
obtained from Eq. p.50p by performing the Fourier trans- 
formation in time variable t ^ t' . We represent it as 
the sum of two terms describing two separate scattering 
mechanisms 



and 



C V ~ sm Im 



e±e 



(3.47) 



The dimensionless parameter has been introduced in 
Eq. (|2.6p . We make an expansion of the integral ker- 
nel 'Kt,t' in Eq. (|3.25p to the first order in CPj^ using the 
relation 



1 -Pf CtT,t2 



Stdis = Stdc + St 



(3.52) 



In p.52p the first term corresponds to the scatter- 
ing off the impurities in the absence of the microwave 
radiation^^ 



Stdc/ 



dip' lyje + W^^,) fje + W^^,,^') ~ /(e, ip) 
2tt I'o T^i-ip 



(3.53) 



(3.48) r^Yie second term in Eq. ([332]) 



-y{e + W^^' ± io) [/(e + W^^. ± v') - /(£, ^)]} (3.54) 
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describes the scattering processes of electrons off impu- 
rities with participation of a microwave radiation quan- 
tum. For the analysis limited to the first order in the mi- 
crowave power, all multi-photon processes are neglected. 
More specifically, the first term in Eq. p.54p represents 
the impurity scattering with one photon emitted (ab- 
sorbed) virtually and can be thought of as the renormal- 
ization of the impurity potential by microwave radiation. 
This term is taken into account in Sec. Ill Al as a linear 
in contribution to jc\ in Eq. (|2.2p . The second term 
in Eq. (|3.54p describes the real processes of emission (ab- 
sorption) of one microwave quantum accompanying the 
impurity scattering. It is this term which produces the 
oscillatory lu/uJc dependence of the magneto-resistance 
and corresponds to ji^ in Eq. (|2.2p . 

We notice that the obtained collision integral vanishes 
in the clean system for any frequency lo away from the 
cyclotron resonance. In the clean limit, the conductivity 
tensor can be found by applying the Kohn's theorem^ 
argumentation. The conductivity tensor of an interacting 
system, which is galilean invariant, is identical to that of 
the non-interacting system. '^^ It follows then that away 
from the cyclotron resonance, uj ^ lOc, the electric field 
appears in the collision integral only in a combination 
with the disorder scattering rate l/rg. 



IV. MAGNETO-OSCILLATIONS IN THE 
PRESENCE OF AC AND DC EXCITATIONS 



In this section we calculate and analyze the dissipative 
current 

j^2evF J ^cosv? J iy{e)f{e,ip)de , (4.1) 

where the distribution function is determined as a solu- 
tion of the kinetic equation Eq. (|3.42p with the collision 
integral given by Eqs. ([X5^ . and ([XM)) . 



A. Solution of the kinetic equation 

We look for the solution of the kinetic equation, 
Eq. (inHSl), in the form 

/(£, ^) = Me) + SUe, ^) + Sfoie) + <5/i(e, ^) , (4.2) 

where the first term is the equilibrium Fermi-Dirac dis- 
tribution function, Eq. (|2.10p . The second term is the 
classical solution corresponding to the constant density 
of states 

Sf ci{£, f ) ^ -defr ^^^"^ cos if , (4.3) 

leading to the Drude result for the longitudinal conduc- 
tivity at large Hall angle, Eq. (|2.13p . 

The third and fourth terms in Eq. (|4.2p are the zeroth 
and first angular harmonics of the correction to the dis- 
tribution function resulting from the quantum oscillatory 



component of the density of states in the collision integral 
(I3.52p . For WcTtr » 1, we keep only the isotropic compo- 
nent and the first angular harmonic of the distribution 
function: 



Sfoie) 
Sfiie.v) 



Xd^fxl sin 



27re 



Al cos 



27re 



(4.4a) 
cos V5 (4.4b) 



The coefficients /, Ai and A2 are fixed by the kinetic 
equation Eq. (|3.42p . The calculation of these coefficients 
is outlined in the Appendix. 

The amplitude / of the isotropic part of the distribu- 
tion function is 



1 



TT T< 



-1 



+ Tg ^ - 7(edc) + 2a'^^7(edc) sin^ Trea 



edc7'(edc) - 27reacSin27reacJ'u;7(edc) 
- 2edc sin^ 7reacJ'i^7'(edc) 



Here, functions 7(edc) and 7(edc) are defined as 

Jni'^'^dc) 



r(edc) = ■ 



and 



7 



;edc) =^ J^(7redc) 



2t, 



n+l 



2t„-1 



(4.5) 



(4.6a) 



(4.6b) 



in terms of the Bessel functions Jn (Tredc) of the order n 
and angular harmonics 1 /t„ of scattering rate off disor- 
der, Eq. (|3.23p . Term Tq"^ in the denominator of Eq. (|4.5p 
is the n = harmonic of scattering rate off disorder and 
coincides with the quantum scattering rate off disorder. 
Coefficient / contains the inelastic relaxation time Tee, 
Eq. p.4ip . Equation (|4.5p was obtained in Ref . [2I for 
arbitrary value of ejc and 7^1 = and in Ref. [iSj for 
Edc 1 and ^ 1 and has been shown to be important 
for the description of MIRO and HIRO at small voltages. 
Equation (14. 5p determines the behavior of the isotropic 
non-equilibrium component of the distribution function 
in constant electric field of an arbitrary strength. 
The amplitude Ai is given by 



Al 



I 



(7'(edc) - 2T„7'(edc) sin^ ne^^) 
+ ^ (edc7"(edc) - 27reac sin 27reacJ'u;7'(edc) (^-T) 
- 2edc sin^ 7reacJ'w7"(edc 



And finally we obtain the amplitude A2 in the following 
form 

^2 - — (7'(edc) - 25'^7'(edc) sin^ Treac) ■ (4.8) 

TTUJr 
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Equations (|4.5p . (|4.7p and (|4.8p determine the first two 
angular harmonics of the distribution function through 
Eq. (|4.4p . This allows us to compute the current density 
as discussed in the next subsection. 



B. Non-linear current 

In this section we calculate the dependence of the cur- 
rent to the first order in power on parameters Cdc and 
Cac, characterizing the strength of the dc and ac excita- 
tions respectively. We substitute the distribution func- 
tion Eq. (14. 2p with factors / and ^1,2 given by Eqs. (|4.5p . 
(|4.7p and (|4.8p to the expression for the dissipative cur- 
rent, Eq. (liTT]) : 



j(edc,eac) =(7uE + Sj^ 



(4.9) 



where the Drude conductivity ctd is given by Eq. (|2.13p 
and the leading correction in to the current has the 
form 



Sj = {Ai - A2) evpyo ■ 



(4.10) 



Substituting Ai,2 from Eqs. pTTj) and p~8| to 
Eq. (j4.10p . we represent the correction to the current in 
terms of the dimensionless function F (edc, fac): 



(5j = 2aD^A2F(edc,eac). 



(4.11) 



The function F(edc, Cac) hi Eq. (|4.1ip describes the lead- 
ing correction to the classical value of the current due to 
oscillations of the electron density of states. This func- 
tion can be represented as a combination of the "dis- 
placement" term ^d(edc,eac) and the "inelastic" term 
(edcjCac), arising due to the non-equilibrium isotropic 
component of the distribution function: 



F (edc7 £ac) — -Pd (fdc, Cac) + {^dci Cac) 



(4.12) 



The "displacement" contribution _fd originates from the 
modification of the scattering rates in crossed electric 
and magnetic fields, while neglecting the effect of electric 
fields on the isotropic component of the electron distri- 
bution function. We have 



Ttr 



7"(edc) 
7r2 

27^ ( 7reac7' sin 27rea 



Edc 



7 Sm TTEac 



(4.13) 



with 7 = 7(edc) and 7 = 7'(edc) defined by Eqs. (|4.6p . 

The remaining "inelastic" contribution is proportional 
to the amplitude / of the isotropic and energy-dependent 
non-equilibrium component of the distribution function 
Sfo{e). This contribution is sensitive to the energy relax- 
ation rate 1/tco at sufficiently small values of edc: 

2 7^ - 27'y^ sin" ^£ac 

Ttr fdc To^ ^ + Tq" ^ - 7 + 2^T^ sin" Treac 

X (cdcl' - 27reac75' uj f^in ^Treac ^e^c"^ J'lj sin '^^ac ) ■ 

(4.14) 



Equations (|4TT|) . (|4T2l) . (|4T3)) and (|4T4)) give the ex- 
plicit expression for the current in response to the ap- 
plied dc electric field with strength E = edc^c/(2|e|i?c) 
in weak microwave fields ^ 1. 

First, we analyze the properties of functions 
^d(edc,eac) and i^i(edc,eac) in a weak dc electric field, 
fdc <C 1. In this case 



7(edc) = 


1 

To 


2rtr 32 T, 


(4.15a) 


7(edc) = 


1 

Ttr 


^^2,2 1 


(4.15b) 


= 1/t„= 


is 


the quantum scattering 


rate off 




1 

Ttr 


1 1 

To Tl 


(4.16) 



disorder. 



is the transport scattering rate written in terms of har- 
monics of scattering rate, see Eqs. (|2.13p and (|3.23p . and 



1 

T* 



3 

To 



4 

Tl 



1 

T2 ' 



(4.17) 



Substituting Eqs. (jil^ to Eqs. (|iT^ and we 

obtain 



Fd = l-^ 

T« 



3 

-TT^e^^ + (-Treac sin 27rea 



sm vreacj 

(4.18) 



and 



Fi = -2 



1 - (Ttr/n)J't^ sin" TTCac 

Ttr/Tco + 7r2e^^/2 -|- 2?^ sin" Trea 



X [7r"e^^ (l - (Ttr/T,)T^ sin" TTCac) + 27rJ'^eacSin27reac] ■ 

(4.19) 



In smooth disorder, Ttr ^ t,. In this case, Eq. (|4.19p 
coincides with the contribution to the electric current, 
considered in Ref. [isl . which we already presented in 
Eq. (|2.20p . The non-linear dependence on the applied 
electric fields occur at edc ^ \/ Ttr/Tco- 

In strong dc electric fields functions 7 (edc) and 7 (edc) 
can be evaluated for edc ^ 1 using the asymptotes of the 
Bessel functions. As a result, the we find the "displace- 
ment" contribution -fd(edc,eac) in the form: 



Fi 4 1 

— =^ sin2^edc(l-2a'„) 

Ttr Tl" edcTTT 

_^ 4 Tc^ y-v edc ± Ca 

± 



edcT^ 



edc 



■ sin27r(edc ± Cac) 



(4.20) 



with l/r^ — cos7m/T„ being the back scattering rate 
off disorder. The function -fd(edc,eac) has an oscilla- 
tory dependence on the parameter edc- We presented 
the corresponding expression for the non-linear contri- 
bution to the current in Sec. |TT1 Eqs. (PTTS)) and (P7T5)l . 
where the same result was obtained from semi-qualitative 
arguments and the stationary phase approximation. 
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For the "inelastic" contribution Fi in the limit edc ^ 1 
we obtain the following estimate 



8 



1 



To 



edcTjr TTredc 



COS 27redc 



1 — ( 4 sin Treac H tan 27redc sin 27rea 



(4.21) 

We notice that this contribution Fi(edc,eac) is smaller 
than Fd(edc, Cac) by factor (1/edc) in strong electric fields 
Edc ^ 1- For a system with smooth disorder -Fi(edc,eac) 
contains also an additional small factor tq/tj^ <^ 1. This 
smallness allows one to neglect the "inelastic" contribu- 
tion to the current in strong electric fields. The first 
line of Eq. (|4.2ip coincides"^- with the asymptote found 
in Ref. 23. The second line of Eq. (|4.2ip describes the 
effect of microwave field in the bilinear response. 



C. Model of disorder 

The purpose of this section is to analyze our results 
for a particular model of disorder. The disorder in 
GaAs/AlGaAs heterostructures is mainly due to remote 
charged donors. The potential created by these donors 
is smooth on the electron wave-length scale. Therefore, 
such potential is characterized by an exponentially sup- 
pressed backscattering amplitude. As we saw in Sec. Ill Al 
the latter is crucial for the onset of non-linear magneto- 
oscillations at higher current densities. On the other 
hand the relatively weak in-plane disorder creates the 
scattering potential of the short range giving rise to a 
finite backscattering amplitude. 

An adequate model should describe a mixed disor- 
der that includes impurities with narrow and wide an- 
gle scattering. ■^^ Specifically, we use the expression for 
the angular harmonics of the scattering rate employed in 
Ref.lii 



1 



Tsh 



1 



1 



1 + xn 



(4.22) 



For this model we have the following expressions for the 
quantum and transport scattering rates 



1 



1 

To 



1 

Tsh 



1 



1 

Ttr 



1 

Tsh 



1 



X 



1 



X 



(4.23) 



We perform calculations under the assumptions: WcTq < 
1 and WcTtr 3> 1, which can be met for disorder described 
by Eq. (|4.22p . provided that x ^ 1 and Tgh ^ Tsm. 
The last inequality means that the long-range disorder 
is stronger than the short-range disorder. The scattering 
time T„ is given by 



1 



Tsh 



(4.24) 



Equations (|4.23p and (|4.24p determine the ratio Ttr/r*, 
which appears in Eqs. (|4.18p and (|4.19p for the non-linear 
contributions to the current at weak electric fields. 

For very small values of edc and T^^, we have the fol- 
lowing estimate 

_ 27ee ^^2^2^ _^ 2TTV^e,, sin2^eac) . (4.25) 

Ttr 

Comparison of this expression with Eq. (|4.18|) shows that 
the linear in and e^^ contributions to the non-linear 
current are dominated by the "inelastic" mechanism if 
Too/Ttr ^ Ttr/T*. If Only loug range disorder is present 
in the system, we obtain for the ratio Ttr/T* — 12x <C 1 
and the "inelastic" mechanism is indeed dominant for 
sufhciently low temperatures, while t^c ^ TtrX ~ tq. 
At higher temperatures the interaction effects suppress 
the oscillations of the density of states and decrease the 
overall amplitude of non-linear current. In the case of 
mixed disorder, parameter Ttr/T* could be of the order 
of unity and the "inelastic" and "displacement" mecha- 
nisms become comparable at lower temperatures, when 
Tee ~ Ttr, but Still Tee ^ Tq and thc osciUations of the 
density of states are not smeared by the interaction ef- 
fects. As a consequence, there exists a regime where the 
linear photo-resistivity may have significantly weak de- 
pendence on electron temperature.'^^ 

We use the model for disorder described by Eq. (|4.22p 
to evaluate the back scattering rate, which characterizes 
non-linear contributions to the current in strong electric 
fields edc S> 1. Performing summation over index n for 
6* = TT in Eq. ^2^, we find 



1 



1 



Tsh 



1 27r 

fsm 



exp 



TT 



(4.26) 



This rate controls the magneto-oscillations for edc ^ 1, 
Eqs. (|4.20p and (|4.2ip . For small values of x, the second 
term can be disregarded. 

To analyze the case of arbitrary edc, we evaluate func- 
tions 7(edc) and 7(edc), Eqs. (|4.6p . in the limit x ^ 1- 
We perform summation over index n in Eq. (|4.6p using 
the Poisson formula and neglect exponentially small cor- 
rections, similar to the second term in Eq. (|4.26p . The 
result reads^^ 



7(edc 



Tsh 



1 



1 



1 



7(edc) =— (Jo(^£dc) " ^i(7redc)) 

Tsh 



(4.27a) 
(4.27b) 



X 



1 



(l + X'T^e^Js/^ 



Equations (|4.27p are convenient to calculate the non- 
linear contribution to the current at arbitrary value of 
edc- The resulting function F(edc,eac), Eq. (|4.12p . is 
plotted in Figs.[Ha,b). We also compare contributions 
^d(edc, Cac) and -Fi(edc, Cac) to the net non-linear current. 
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Both functions are shown for the same set of parame- 
ters in Fig. [5]^c). The "inelastic" contribution _Fi(edcj Cac) 
has a large variation at small ecic7 and vanishes at larger 
values of edc- On the other hand, the displacement com- 
ponent -Fd(edcjeac) is an oscillatory function of the pa- 
rameter edc with only weakly decaying amplitude. 




FIG. 2: (Color online) The function F(e(jc,eac) as defined by 
Eqs. (|4.12|l . (|4.13|) . and (|4.14|) with 7(edc) and 7(edc) gi ven by 
Eqs. (I4.27P for the model of disorder introduced in Sec. IIV Cl 

(a) r,„, = 1, r,h = 1, r^o = 0.5, x = 0.0001 and 7^ = 0.25 ; 

(b) Tsm = 0.1, Tsh = 1, Tec = 10, X = 0.001 and = 0.01. In 
both cases the thick solid (red) line, the dashed (blue) line, the 
dashed-dotted (black) fine and the thin solid (black) line are 
used for eac = 2.75, eac = 3, eac ~ 3.25, and eac ~ 3.5, respec- 
tively. The panel (c) shows the inelastic contribution, Fno, 
solid (red) line, and displacement contribution F^, dashed 
(blue) line for the set of parameters of panel (b) at eac ~ 2.75. 



For edc ^ 1, when the "inelastic" contribution becomes 
small, the function i^(edc,eac) coincides with the "dis- 
placement" contribution Fd(edc, Cac)- The latter is a sum 
of two terms. One term contains 1/Tsm, this term varies 
smoothly as a function of edc on the scale Cdc ~ 1/VX- 
Another term contains l/rgh and oscillates as a function 
of edc with period equal to unity. The corresponding os- 
cillating term can be written as 



Ttr 



Edc 



(Jq (Tredc) - Ji (Tredc))" sin^ 7re 
(Jo(7redc) - Ji(7redc))'sin27re, 



(4.28) 



This equation indicates clearly the vital importance 
of the short-range disorder, characterized by the finite 
backscattering rate l/rgh for the onset of the magneto- 
oscillations in the non-linear transport regime. 



D. Differential magneto-resistance 



Now we apply Eq. (|4.9p for the electric current response 
to the applied dc electric field to describe the longitudinal 
differential magneto-resistance 



(4.29) 



In Eq. (|4.29p we have chosen to be the current direc- 
tion, so that is the electric field component parallel 
to the current. In the limit lUc S> I/th- we write 



(4.30) 



where the electric current density j is given by Eq. (|4.9p 



PD] 



(4.31) 



calculated with a total electric field E ~ paj, where pn 
ujc/{e'^Vpi'o) is the Hall resistance and 



2|e|(pHj)i?c 



■^dc 



47rj 
epFUJc ' 



(4.32) 



We finally obtain for the oscillatory part of the differen- 
tial resistance, Sp = p — pjj the following expression 



PD 



"^dc 



^dc^(^dc' ^ac) 



(4.33) 



with function F(e^^, eac) given by Eqs. ([i7T^ . ([iTT^ and 
(|4.14[) . Figure [3] shows the differential resistance found 
for the specific disorder model considered in Sec. IIV Cl 

We analyze our results in the limiting cases of small 
and large dc currents. We start with the discussion of the 
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FIG. 3: (Color online) The upper and lower panels show the 
differential resistivity as a function of ejc for values of param- 
eters used in Figs.[2|a) and (b), respectively. The differential 
resistance is obtained by substitution of F{e^^, eac), shown in 
Figs.IHa) and (b), to Eq. (|433|) . 



differential resistance Sp at small direct current, ^ 1' 
when in samples with smooth disorder the non-linear be- 
havior originates from -F](e;^^, eac) contribution. Keeping 
only this contribution in Eq. (|4.33p . at e^^ = we have 



PD 



= - 2y 



27r£ac sin 27reac 

Ttr/Too + Sin^ TTEa 



(4.34) 



As e^^ increases, the differential resistance decreases from 
the above value on the scale e'^^ 
a non-monotonic behavior at e^^ > 1. However for large 
values of e^^, the "inelastic" contribution to the differen- 
tial resistance vanishes fast, as discussed in Sec. lIVBl 

We now analyze the differential resistance at large val- 
ues of direct current, e^^, ^ 1. In this limit the "dis- 
placement" contribution dominates over the "inelastic" 
contribution. Therefore, we keep only -^^(edc' ^ac) and 
find 



of Ref.l23lat ?^ = 0. At finite 7^ the factor (1-25'^) rep- 
resents the impurity potential renormalization due to the 
virtual absorption and emission of photons. These radia- 
tive corrections tend to suppress the dark resistivity and 
can be interpreted as motional narrowing or averaging 
out the electrostatic potential of impurities in the pres- 
ence of an oscillating electric field. The other two terms 
in Eq. (|4.35p describe the effect of the combined scatter- 
ing off disorder in mixed constant and oscillating electric 
fields when the real microwave photons are absorbed or 
emitted. 




-ac 



-ac 



FIG. 4: (Color online) The differential resistance obtained 
from Eq. (|4.35|l for = 0.25 presented on a grey-scale plot. 
The bright areas correspond to a higher value of the resis- 
tance, (left) Away from the main diagonal in the (eac,e^j.) 
plane the differential resistance has minima and maxima de- 
scribed by the conditions (|4.36p . (right) Close to the main 
diagonal eac — e^^, the differential resistance becomes a func- 
tion of the sum eac + e^^. 



We briefly discuss the properties of Eq. ()4.35p as a 
function of two parameters e^^ and eac- It is plotted in 
Fig. [4] as a grey-scale contour plot. This function exhibits 
a series of maxima and minima in (eacCdc) plane. In 
general, for eac > Ejc local maxima and minima of 
the function defined by Eq. (|4.35p are located at 

l,tac, e^jj,; 



■^dc 
J \ min 



ac , Cjjj, 



(m± l/4,n=F 1/4) 
(m± l/4,n± 1/4) 



(4.36) 



\/ Ttr/Tee and exhibits We further notice that in the region of parameters 



^ Cac + Cdc ' namely not to far from the main 
diagonal in the two-dimensional plain (cac, Cdc) the result 
(|4.35p reduces to 

Sp ^_^ (4A)Vtr 
PD TTT^ 



Sp_ ^ (4A)Vtr 
PD TTTtt 



(1 - 2y^) cos 27re^^ + 2?^ cos 27r(e; 



3 

dc ■ 



(l-2J'^)cos27redc 



(4.37) 



-I- 23'(j cos 27redc cos 27reac sin 27redc sin 27rea 

V Cdc 

(4.35) 

The first term in square brackets coincides with the result 



Although a direct addition of two different parameters e^^. 
and eac has no physical meaning, we believe that, in the 
considered region edc — Cacj it is the second term which 
is responsible for apparent structure^i of the differential 



resistance as a function of the sum eg 



""dc 



see Fig.|3Kb). 
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V. NON-LINEAR MAGNETO-RESISTANCE 
BEYOND THE BILINEAR RESPONSE IN 
MICROWAVE FIELD 

In this section we extend the analysis of the magneto- 
transport at edc ^ 1 to the case of arbitrary strength 
of microwave radiation. We again consider the case of 
relatively high microwave frequency uj 3> l/Ttr- The last 
condition ensures that the oscillating in the variable {t + 
t')/2 part of the distribution function is small compared 
to the stationary one. We therefore can assume that 
the distribution function is homogeneous in time, ft^v = 
ft-t'- Beyond the bilinear response, Eq. (|3.50p . obtained 
as the expansion of Eq. (|3.25p in powers of T^, is no 
longer valid. We perform the time averaging of the exact 
collision kernel Eq. p.25p over the microwave oscillation 
period: 



dip' 

"2^ 



(5.1) 



X e 



jpF(n,^-ii^/)Cf't, 



For the circular polarization the displacement due to the 
microwave field Ql, is given by Eq. ()3.47|) . The time 



averaging in Eq. (j5.1[) results in 



with Jq{x) being the Bessel function and 



sin 



(5.2) 



(5.3) 



The time averaged collision kernel takes the form 



(5.4) 



X Jo I^Q,^-,^' sin ■ ^ 

The substitution of Eq. ((O)) to Eq. ((33T|) leads to the 
following expressions for the "out" collision term 

Stout{/}t,t' = — / -^-^ — 



X \e 



dip' A . 

-z Jo (Qip~ip' siuTreac) 

W iT,. r , -iW iT,. r 



The "scattering-in" term Eq. (|3.36[) reads 

J l-K T^-^i 

X -'^o I Qip-p' sm 

X (/t,t' - ^ft-T^,t' — ^ft.t'-T^ 



(5.5) 



(5.6) 



We rewrite the kinetic equation Eq. (|3.42p for homo- 
geneous in time distribution function in the energy rep- 
resentation 



d 



(5.7) 



The full collision integral for electron scattering off dis- 
order, Stdis{/e;i^}e;i,3, is a sum of the scattering "in" and 
"out" terms, given by Eqs. ((53| and ((O)) . see Eq. ((3?28)l . 

We discuss the non-linear transport regime ejc 3> 1. 
In this case it is sufficient to consider the "displace- 
ment" contribution to the non- linear current. To cal- 
culate the "displacement" contribution to the electric 
current, we substitute the equilibrium distribution func- 
tion Eq. (|2.10p to the collision integral in Eq. (|5.7p . Com- 
bining Eqs. (|5.5p and (j5.6p and writing the result in en- 
ergy representation, we obtain 



Stdis{/T(£)}s;v 



The first kernel, 



dip' Ko ~ 2ARe{e*('^+^^^')^'^ifi} 



27r 



(5.8) 



Ko = Jo [Q^^^' sin -^j Me + W^^,) - Me), (5.9) 

is the part of the collision integral to the zero order in A. 
In particular, it describes the classical Drude conductiv- 
ity. Although the microwave radiation complicates the 
form of Kq, we notice that this kernel can be represented 
in terms of the series expansion in powers of d/de, ap- 
plied to the Fermi distribution function. The observables, 
such as electric current, are determined by energy inte- 
grals. Because the integrals of all derivatives of frie) 
of the second order and higher vanish, we do not expect 
any effect of microwave radiation within our model on 
the conductivity to the zeroth order in A. 
The second term in Eq. (|5.7p is 



Ki ^Jq i Q^-^' sin 



— TTCac 



- Jo {Q^-v' siuTreac) frif-)- 



)Me + W,^^,) 

(5.10) 



In the high temperature limit, T 3> \e\ERc and T ^ w, 
Eq. (|5.10p can be reduced to 



dMe) 



Jq iQf-f' siuTreac) 



Ji [Qip-if' smTreac) cosTrea 



(5.11) 



Similarly to the case of the small power, we look for the 
correction to the distribution function in the form 

2tt£ 

6f{e, (fi) = 5fc\{e, (p) + XdefrM cos cos ip . (5.12) 

Substituting 5f{e, ip) to Eq. (|5.7p . we find that (5/ci(e, ip) 
is given by a sum of Eq. (j4.3p and higher order derivatives 
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of fris), as discussed below Eq. (|5.9p . The term of the 
first order in A determines the value of coefficient Ai: 




-0.2' ■ ■ ' 

3 4 5 6 



FIG. 5: (Color online) The function Ai evaluated by (a) nu- 
merical integration in Eq. (|5.13p . solid line (blue), and (b) 
by the stationary phase approximation applied to Eq. (|5.13|l , 
dashed line (red) for = 8.25, eac = 2.25. 



Our analysis in this section is applicable for ejc ^ 1, 
when we can utilize the saddle point approximation, to 
perform angular integrations in Eq. (|5.13|) . Angles ip and 
(/?' which make the phase of the integrand stationary are 
if = ±7r/2, if' = ^7^/2 and correspond to 2Rc jumps 
along the electric field direction, see Sec. |TT1 For ^/TZ ^ 
Edc the angular integrations in Eq. (|5.13p results in 



Ai = ^ — sin27redc-^o ( 4v J'w sinTreac 

^ — cos 27redc cos TTEac v^5\^ Ji f4\/T^ sinTreac") • 

TT^T^ edc ^ ^ 

(5.14) 



The result of numerical evaluation of the coefficient Ai, 
Eq. (|5.13p . and its approximation, Eq. (|5.14p . arc shown 
in Fig. [5] The stationary phase approximation progres- 
sively improves as the parameter edc increases. 

The non-linear contribution to the electric current is 
given by Eq. (I4.10p in term of coefficients Ai, Eq. (|5.13p . 
and A2 = 0. Representing the non-linear contribution to 
the current in the form of Eq. (|4.1ip . we write 



i^(edc,eac) -Ai. (5.15) 

edc 

The function F(edc,eac) found from Eqs. (|5.14p and 
(|5.15p reads 



4rt, 



sin27redc^o (^4V5\^sin vrCac^ + cos 27redc cos Treac V^P^ Ji (i^/V^ sinnea 



(5.16) 



The oscillatory correction to the differential magneto-resistance is obtained by substituting Eq. (|5.16p to Eq. (|4.33p : 



SpU) _(4A)Vt, 



PD 



COS 



27redc^o 



Ax/y^ sin? 



— — sin 27redc cos TrCac \P^J\ ( 4-\/ sin Trea 

Edc 



(5.17) 



The relation (|5.17p reduces to Eq. (|4.35p obtained earlier in the weak power limit as expected. 



VI. DISCUSSION AND CONCLUSIONS 



In this paper we have presented a comprehensive the- 
ory of an out of equilibrium two-dimensional electron sys- 
tem (2DES) in a magnetic field in the case when the dis- 
ordered potential has finite scattering amplitude on an 
arbitrary angle. Then, we applied this theory to ana- 
lyze electron transport in the presence of constant and 
oscillating in-plane electric fields. We showed that the 



electric current has an oscillating component as a func- 
tion of the strength of the constant electric field and of 
the frequency of an oscillating electric field. We have 
investigated the position of maxima and minima of the 
differential resistance on the plane of parameters Cdc and 
eac and found a qualitative agreement with experiments 



ofRefs. [mi]. 



Actual value of the differential resistance as a func- 
tion of the current through a 2DES sample depends on 
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a number of different scattering rates, including the full 
quantum scattering rate, the transport scattering rate 
and the backscattering rate. Our analytical results may 
be applied to the experimentally measured differential 
resistance to evaluate these scattering rates and to ob- 
tain a detailed picture of the origin and structure of the 
disorder in high-mobility 2DES samples. In particular, 
we would like to emphasize that the behavior of the dif- 
ferential resistance at large currents in experiments of 
Refs. [T9ll2^l25j suggests that the disordered potential 
has a noticeably strong short-range component, respon- 
sible for the finite backscattering rate. 

The presence of the short range disorder also modifies 
the current behavior in weak constant electric fields and 
at weak power of an oscillating electric field. In partic- 
ular, there are two competing contributions to the cur- 
rent. One contribution originates from the modification 
of the isotropic component of the electron distribution 
function by electric fields, we refer to this contribu- 
tion as an "inelastic" contribution. The other contri- 
bution is due to the modification of electron scattering 
rate off disorder by electric fields, known as the displace- 
ment contributioniiSiii^ii^iS While in smooth disorder in 
weak fields the "inelastic" contribution dominates in the 
whole range of temperatures at which the non-linear cur- 
rent is expected to survive, the short range disorder may 
make these two contributions comparable. As a result, 
the strong temperature dependence of the non-linear cur- 
rent in weak fields is expected only in samples with suf- 
ficiently weak short-range disorder. 

We also studied the dependence of non-linear current 
on the applied power of oscillating electric field. The 



expression (|5.17p obtained in the non-linear transport 
regime ejc S> 1 is not limited to the small microwave 
radiation power. This expression shows that in strong 
microwave fields the non-linear current has a different 
dependence on parameters eac and Cdc, Eq- (|1.1|) . which 
qualitatively differs from the corresponding expressions 
in the weak power limit. We believe that the study of 
non-linear current at strong radiation power may bring 
additional opportunities for study microscopic character- 
istics of high mobility electron systems. 
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APPENDIX A: CALCULATIONS OF THE 
ELECTRON DISTRIBUTION FUNCTION 

In this appendix we give the details of the calculation of 
the coefficients /, Ai and A2 of the distribution function 
in the form of Eq. (|4.4p . Substituting Eq. (|4.4|) to the 
kinetic equation (|3.42p with the collision terms specified 
by Eqs. (|3.53p and (|3.54p we obtain the following system 
of equations 



1 



1 



Toc(e) To 



sin^^ :^/K^^f{e) - 2'^^ 'Kta.cK^p^p,{e) + LJsin27reacM;^;p/(e) 
\ 

/ {M^^,{e) - sin^ Tre^^M^^,{e)] ) , 



(Al) 



and 



— - — cos =( sinip\K^^'[e) — 27'^ [sin neacK^^iye) — a;sin27reacM";p;^'(e)J 

+ I [M^^,{s) - 2T>^ sin^ Tit^^M^^.{e)\ } 



=/( sin cos 



sme 



2^(£ + W^^. ± ^) , 
cos ^ ( M^^> (e ± w) 



sine 



± L 



cos 



r 



(A2) 



(A3) 



where (. . .) stands for the averaging over angular vari- 
ables if and ip' . Here we have introduced integral kernels 



The kernels K^p^pi{e) and M^p^pl{£) are obtained from 
Eq. (IA4p by replacing r ^ f as given by Eq. (|2.5p . 



K^^p, (e) = - 2 cos 



2TT{e + W^^,) W^p. 



Tip — ip' 

sin [2n{e + W^p^pl)|w^ 



(A4) 
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Keeping only the energy independent part in Eq. (|A3p 
surviving the subsequent energy integration we rewrite 
it as 



/ (sin if [M^^, (0) - sin^ ne^^M^^, (0)] ) . 

(A5) 



For the angular averages appearing in Eqs. IjAip . (IA2p . 
and (|A5p , we have the following expressions 

/ \\ 2eERc I, , . 27re , ^ r. \ 

(K^^i (e)) = 7 (edc) sm , (A6a) 



27re 

{M^^, (e)) = 7(edc) sin , (A6b) 



(sin ipK^p^pi{e)) 



5 cos 



, (A6c) 



{sm(pMu:,ip'[e)) — — cos 



— . (A6d) 



In Eqs. (|A6p the function 7(edc) has been introduced in 
Eq. (|4.6ap . The angular averages involving bared func- 
tions K and M are given by Eqs. (jA6P with 7(edc), 
Eq. (|4.6bp . replacing 7 (edc)- Relations (jA6p allow us to 
perform the angular integrations in Eqs. (jAip . (|A2p and 
(|X3l) leading to the expressions (Ii3)l , ((i?7)) and (g^ of 
the main text. 
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